Multifractality of quantum wave packets 
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We study a version of the mathematical Ruijsenaars-Schneider model, and reinterpret it phys- 
ically in order to describe the spreading with time of quantum wave packets in a system where 
multifractality can be tuned by varying a parameter. We compare different methods to measure the 
multifractality of wave packets, and identify the best one. We find the multifractality to decrease 
with time until it reaches an asymptotic limit, different from the mulifractality of eigenvectors, but 
related to it, as is the rate of the decrease. Our results are relevant to experimental situations. 

PACS numbers: 05.45.Df, 05.45.Mt, 71.30.+h, 05.40.-a 



Multifractal properties have been characterized in sev- 
eral physical contexts, from turbulence [l| to the stock 
market or cloud images Similar features were 
also recently observed in quantum mechanics, or more 
generally in complex wave systems, in specific regimes. 
In particular, multifractal wave functions can be seen in 
condensed matter systems, where they are believed to 
characterize electrons at the Anderson metal-insulator 
transition [Jill- These properties were also observed in 
quantum Hall transitions [7|, in Random Matrix models 
[1, [§], and in other physical models l3"12|- They are 
also visible in the different context of pseudointegrable 
systems, for which constants of motion exist, but dy- 
namics takes place on manifolds more complicated than 
the tori characteristic of integrability [13|-|20| . 

Many theoretical studies have been devoted to these 
quantum multifractal systems, and several important 
results have been obtained. In parallel, experimental 
progress opens the way to direct observation of multi- 
fractality; hints of such properties were seen in waves in 
elastic media 21 1, disordered conductors [l^] and cold 
atoms [i^. The theoretical studies were mainly concen- 
trated on eigenvectors (stationary states) of the physical 
systems considered [4l420l|. In contrast, the experimen- 
tal protocols in general involve the propagation of Wave 
Packets (WP), for which results on eigenvectors are a 
priori not directly applicable. In order to further charac- 
terize experimental results and interpret them, it is there- 
fore important to have detailed results on the multifractal 
properties of WP of such systems, and to understand how 
they relate to the multifractality of eigenvectors, which 
was the focus of most of the theoretical results. There 
have been some attempts in the past to relate global 
properties of WP, such as spreading laws or shape of 
envelopes, to the multifractal properties of eigenvectors 



or eigenspectra [2J-|28[. Other works found specific ex- 
amples of multifractal WP [2^,[3^. However, the general 
existence and origin of the multifractality of WP is still 
unclear. In this Letter, we use a mathematical model 
known as the Ruijsenaars-Schneider model [sij and rein- 
terpret it as the quantization of a pscudo-integrable map 
on a phase space whose size can be varied. The prop- 
erties of the eigenfunctions can be continuously tuned 
through system parameters from a weak to a strong mul- 
tifractality regime, enabling us to systematically compare 
multifractality of WP and eigenvectors in the full range 
of possible regimes. Our computations show that several 
methods can give different results for measuring this mul- 
tifractality of WP, even though they are equivalent for 
eigenfunctions, and they allow us to identify the optimal 
method. Our numerical and analytical results show that 
one can systematically relate multifractality of WP and 
its time evolution to the one of eigenfunctions, opening 
the possibility to probe these properties in detail through 
experimental observations. We present here the main re- 
sults, details will be published elsewhere. 

We consider a periodically kicked system with period 
T and Hamiltonian H{p, q) = ^ + V{q) Y,n ^i* " ^T), 
with potential V{q) = — 7{?} (here {q} denotes the frac- 
tional part of g), and {p,q) are the conjugated momen- 
tum and position variables. The classical motion takes 
place on a cylinder, with the classical map over one pe- 
riod p = p + 7; Q — 1 + Tp (mod 1). The quan- 
tization of this map yields a unitary evolution operator 
given hyU ^ c^'^p /i2h)^-iv{q)/h . in the choice 

of parameters and of the effective constant 2TTh = N 
led to a quantum map on a toroidal phase space inde- 
pendent on the dimension of the Hilbert space N. In 
order to allow for long spreading times for a WP, here we 
choose h fixed and truncate the phase space by taking 
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p € [0, 2TTNh[, or equivalently integer indices P defined 
by p = 2nPh such that 0<P<7V-1. We thus define 
a quantum map over a phase space whose classical size 
grows with N. The evolution operator then becomes U = 
^-i^TP'p-i^2^nQ/Np ^ ^j^i^ ^ e^piiPQ/N)/y/N, 

which yields 
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TV 1 - g2j7r(P-P'+7)/Af 



(1) 



This system corresponds to the Ruijsenaars-Schneider 
map, well-known in mathematics since it describes the 
Lax matrix of certain many-body classical systems [l9l . 
[20| . This model reinterpreted physically in this way has 
many advantages. The multifractality of eigenvectors has 
been shown to depend on the parameter 7 in a continu- 
ous way, enabling to probe all the regimes from weak to 
strong multifractality [isl - fioj . In addition, the simplicity 
of this ID model makes analytical calculations and nu- 
merical computations tractable. In the numerical results 
below, we choose to replace the kinetic term in ([T]) by 
random phases <f>p, in order to get averaged quantities 
while keeping the same physics. 

We consider the evolution of a WP initially localized on 
one single momentum state, ^fp-* = S{P — Pq). Iterations 
of the map make the WP spread out. Analytical calcula- 
tions are possible in the regime of small t and 7 = /c + e 
where k is the nearest integer and e is small. Using a 
tailored version of perturbation theory, originally devel- 
oped for eigenvectors we were able to compute the 
average WP over random phases <I>p. Long and techni- 
cal calculations which we do not present here yield a very 
simple general formula in this regime for the average WP 
for P 7^ Po - kt: 
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It implies that close to integer values of 7, the WP dis- 
plays a single peak moving at speed k. Actually, formula 
([2|) is close to the numerical results even for quite large 
values of t and for 7 far from integers, provided k is re- 
placed by 7 and e by sin(7r7)/7r. This can be seen on the 
insets of FiglTJ where the average wave function is shown 
together with the formula ([2]) for three different values of 
7. Discrepancies can be visible only by zooming close to 
the center of the distribution. In the case 7 ~ 0.5, one 
can actually distinguish two peaks, one staying at the ini- 
tial position, the other one moving faster than 7, but the 
tails are still perfectly reproduced by Eq.(I2]). The peak 
at the origin for small 7 can be interpreted as a manifes- 
tation of strong multifractality, since it can be related to 
the correlation dimension 2J, |25|, |27 1 . 

Several regimes can be characterized in the evolution 
of the WP. At t = 0, the WP is locahzed at P = Pq. As 
time increases, the WP spreads according to Eq.(I2|) until 



it reaches the system size. We have checked that the lim- 
iting result for very long times (t = 00 limit) is equivalent 
to the one obtained by diagonalizing the evolution oper- 
ator, replacing the eigenphases by random numbers, and 
transforming back to the momentum basis. The speed at 
which the properties converge to this asymptotic regime 
depends on the multifractality (see below). 

Different methods can be used in order to compute 
multifractal exponents. They are all equivalent for math- 
ematically defined multifractal measures. However, our 
system is discrete and therefore cannot reach arbitrar- 
ily small scales. Besides, WP can be described as a 
smooth average Eq.Q, with superimposed fluctuations. 
Our study focuses on the multifractality of these fluctu- 
ations. However, the presence of a nontrivial envelope 
can affect the computation of the exponents. This ex- 
plains why in our case different methods may give differ- 
ent answers. We thus tested four different algorithms for 
computing multifractal exponents. The first one (mo- 
ment method), widely used in a quantum context, e.g. 
in 3: 01 1 consists in computing the moments of the wave 
function Vq ~ J2p=i I^Pp' for different system sizes TV; 
the multifractal exponent Dq can be obtained from the 
scaling of these moments with through (Vq) oc N^'^i 
with Tq = Dq{q—1). However, this method assumes scale 
invariance of the system as A'^ increases. In the system we 
study the envelope has a nontrivial scaling with N, and 
this effect is hard to disentangle from the multifractal- 
ity due to the fluctuations. Thus while in most regimes 
we found this method to give results equivalent to other 
ones, in some cases it gives nontrivial multifractal expo- 
nents even for a smooth WP (e.g. for the average WP of 
Eq.@). 

We therefore investigated alternative methods, which 
use only one system size N. One possibility is to eval- 
uate the scaling with box size of the moments Vq°^ = 
I]boxos(I]pGbox I^Pp)'* of the wave function summed up 
inside boxes of same sizes (box-counting (BC) method) 

, . Another method uses the scaling of the sum of 
the local maxima of the wavelet transform of the wave 
function at each scale (wavelet method) (32- 34 1. Finally, 
we investigated a method suited for WP spreading with 
time, similar to the moment method but using different 
times instead of different system sizes (time method) [s^] . 
In this case, we study the scaling of the moments of the 
wave function as a function of time. 

Our computations showed that for strong multifrac- 
tality BC and wavelet methods gave the same results. 
However, in the weak multifractality regime, where the 
scaling is strongly dependent on the scales at which we 
choose to fit the data, the BC method appeared to be 
more reliable. We checked that if the envelope ([2]) was 
eliminated, keeping only the fluctuations of the WP, the 
results of the BC method remained the same (data not 
shown). We therefore choose throughout the paper the 
BC method for the numerical computation of the expo- 
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nents. An average was made over the position of the 
center of the boxes, to ehminate a threshold effect hnked 
to the relative position of the WP and the boxes. The 
time method was also explored as it is somewhat natural 
in the context of the spreading of the WP, but in our case 
the method is cumbersome since it requires some knowl- 
edge of the spreading law to be implemented, and less 
reliable than the BC method. For relatively large values 
of q {q > 1) results were similar to those obtained by the 
other methods. 

We recall that the exponents Dq arc positive and de- 
crease for g > from Dq = 1; a.t a fixed value of q > 
the smaller Dq is the stronger the multifractality is. For 
an ergodic wave function one has Dq = 1 for all q. In 
systems such as ours where an average is made over wave 
functions one can distinguish two sets of multifractal ex- 
ponents 0, 0|. For the BC method, the first one corre- 
sponds to {Vq°^) as defined above, and yields exponents 



D„, while the other one corresponds to {\n'P}'°^ 



and 



q, .....^^ v^x.^ , q 

yields exponents D^^^. In cases where moments arc dis- 
tributed according to power laws with small exponents, 
the two quantities can be different, the latter being the 
typical value of the moment for the bulk of the wave 
functions considered, while the former could be domi- 
nated by rare wave functions with much larger moments. 
As the quantity Dq is the most accessible to analytical 
methods, and the most widely studied in the literature, 
we concentrate on it. We have nevertheless checked that 
our results are similar for both quantities in the regimes 
we considered. 

We now turn to the discussion of multifractal expo- 
nents Dq for WP in different regimes of 7. Studies of 
eigenvectors of the map ([T]) have shown 1^ IS] that their 
multifractality is the strongest close to 7 = and the 
smallest for 7 close to nonzero integers. In Fig (T] we show 
the exponents Dq for WP at i = 100 and t = 00 and for 
eigenvectors in different regimes. As long as the WP re- 
mains localized, the multifractal exponents are extracted 
from scales smaller than the typical WP size. Indeed, 
above this scale, the weight is concentrated only in one 
box and docs not depend anymore on the box size (see 
insets of FiglT] where a saturation of the moments can be 
observed). For very small times there are very few scales 
from which to extract the exponents, which may affect 
the precision. 

The comparison of the different curves in FiglT] (top) 
shows that the regime of small 7 corresponds to a strong 
multifractality of both WP and eigenvectors. In this 
regime, it is possible to use a specific pcrturbativc ap- 
proach in order to obtain multifractal properties of eigen- 
vectors [23|. It predicts that for g > 1/2 the multifractal 
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while for q < 1/2 it is 
The results displayed in 



1 ~ q-l ' ^ ' ^{q-l)r{~q)- 

FiglT] (top) show that this formula is also quite close to 
the multifractal exponents of WP even for values of 7 



Da 



Da 




FIG. 1: (Color online) Multifractal exponents Dq vs q for 
7 = 0.05 (top), 7 = 0.5 (middle) and 7 = 0.95 (bottom). 
Red solid curves: WP at f = 100; blue dotted curves: t = 00 
limit; green dashed curves: eigenvectors of (|T]); black dotted- 
dashed curves: analytical theory for eigenvectors (see text). 
Right insets: average WP (initially localized at P = 0) (dotted 
yellow curve) and analytical formula ^ (black solid curve). 
The Dq of WP have been extracted from 1000 random phase 
realizations of size A'^ = 2^® by the BC method applied on the 
2^^ central components with boxsizes from 1 to 64. The Dq 
in the t = 00 limit and for eigenvectors have been extracted 
from 2^'^ vectors of size 2^^ by the BC method with boxsizes 
from 8 to 512. Shaded areas indicate standard error in the 
least-square fit. Left insets: examples of fit (BC method). 



as high as 0.05. Our explanation is that the eigenvec- 
tors and initial WP being both very localized, only a few 
eigenvectors contribute to the WP. Thus in this regime 
the multifractal exponents of the WP yield a direct infor- 
mation on those of eigenvectors. Wc also note that the 
t = 00 limit is reached already for t = 100. 

When 7 goes further away from zero, the multifractal- 
ity of eigenvectors decreases. As can be seen from the 
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numerical data displayed in Fig[T] (middle) for 7 = 0.5, 
the multifractality of WP also becomes weaker. In this 
regime the multifractality for t = 100 and t = 00 are quite 
close, showing that the t = 00 limit is reached quite fast. 
This is all the more remarkable since as shown in the in- 
sets of FigUtop and middle, the WP a.t t = 100 remains 
on average quite localized. Thus even though the enve- 
lope at f = 00 is flat (data not shown), the fluctuations 
which correspond to the multifractal properties are sim- 
ilar. In this regime, the asymptotic limit is thus quickly 
reached, and corresponds to a multifractality weaker than 
for eigenvectors. Our interpretation is that as eigenvec- 
tors are more delocalizcd than for 7 ~ 0, the initial WP 
has significant components on more eigenvectors, which 
lead to an overall decrease of the multifractality as time 
evolution mixes these eigenfunctions. 

When 7 increases and gets close to nonzero integer k, 
eigenvectors display weak multifractality. This can be de- 
rived analytically since the perturbative approach yields 
the expression = 1 — (7— k)'^q/k'^ In this regime 
the WP at t = 00 are also weakly multifractal, as can 
be seen from the numerical data displayed in Fig[T] (bot- 
tom). The two curves are quite close but in this regime 
of very weak multifractality eigenvectors are slightly less 
multifractal than the t ~ 00 limit. For t = 100 the mul- 
tifractality of WP is quite different from the asymptotic 
one at t = 00. We will come back to this point below. 
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FIG. 2: (color online). Density plot of Di as a function of 7 
and time t. Colors denote multifractality strength from white 
(strong) to red/gray (weak). The exponents Di have been 
extracted from 100 random phase realizations of size A'^ — 2^^ 
by the BC method with boxsizes from 1 to 64. 

The global picture for the multifractality of WP is sum- 
marized in FiglJl which displays Di as a function of time 
and 7. The three regimes can be clearly distinguished, 
both in the average multifractality and in the speed with 
which the asymptotic regime is reached. In order to shed 
more light on the way multifractality evolves with time 
we show in FiglS] the time evolution of Di and D2 for 
three different values of 7. While the asymptotic regime 
is reached very quickly for 7 w 0, the rate of convergence 
decreases with 7, as can be checked more quantitatively 
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FIG. 3: (color online) Multifractal exponents Di (top) and D2 
(bottom) vs t, extracted from 1000 random phase realizations 
of size A'^ = 2^^ by the BC method with boxsizes from 1 to 
64. Red/bottom curves: 7 = 0.05; blue/top curves:7 = 0.5; 
green/middle curves: 7 = 0.95. Upper inset : time r (defined 
by -Di(t) — _Di,as/2) vs 7. Here -Di.as is to the mean value 
of Di (t) in the time interval [3000, 5000] . The solid curve is 
the best fit of the form exp {A/{B - 7)) with A = 0.749 and 
B = 1.103. Lower inset: mean second moment vs time for 
7 = 0.95 (green squares), 0.05 (orange crosses), 0.99 (magenta 
circles) and 0.01 (blue triangles). Black dashed lines is the 
analytical formula for the second moment (see text). 



with the data shown in the top inset. Our interpretation 
of this phenomenon is the following; as 7 increases, we 
have already noted that the initial WP has significant 
components on more and more eigenvectors. One can re- 
late the time at which the asymptotic regime is reached 
to the inverse of the average spacing between eigenphases 
corresponding to these eigenfunctions which contribute. 
For 7 « the eigenphases are close to the random vari- 
ables <I>p and only few eigenvectors contribute, leading to 
a very large mean spacing and thus a very short conver- 
gence time. As 7 increases, this mean spacing decreases 
and the convergence time increases. 

A perturbative method similar to the one used to ob- 
tain can be developed for the second moment of 
WP for 7 close to integers. It leads to (Epl*^^!^)"! w 
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1 + 27r2e2(iV2 _ i)/(3iv2)i, valid for t < 3/{2Tr^e^); data 
in Fig|3] (lower inset) confirm the increasing range of va- 
lidity of this formula when 7 gets closer to integers. It 
predicts D2 ~ 0, compatible with the asymptotic limit 
for 7 « 0, but not for 7 close to other integers. This 
analysis therefore further confirms that as 7 gets closer 
to nonzero integers (weak multifractality regime), the 
asymptotic behavior should take longer and longer times 
to appear. 

We have studied the multifractality of individual WP 
in a periodically kicked system through a combination 
of numerical and analytical works. We have compared 
different methods to define and measure it, and assessed 
their usefulness in this context. The multifractality of 
WP was shown to typically decrease with time until it 
reaches an asymptotic limit, which corresponds to the 
model with randomized eigenvalues. This asymptotic 
multifractality is different from the one of eigenvectors 
more commonly studied, but is related to it. The rate 
at which the asymptotic limit is reached can also be re- 
lated to the multifractality of eigenvectors. Our results 
should be applicable to other models and to the analysis 
of experimental situations. 
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